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ON ■ Abstract 
O 

We study perturbative QCD corrections to deeply virtual Compton scattering on an unpolar- 
ized nucleon target in the flavor non-singlet sector to next-to-next-to-leading order accuracy, 



in 
o , 

, restricting ourselves to the kinematically dominant amplitude. The difference between the stan- 
dard minimal subtraction and the conformal scheme, in which conformal symmetry is manifest, 
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(D \ is studied to next-to-leading order. Beyond this order we employ conformal symmetry for the 



evaluation of perturbative corrections. Within a certain class of generalized parton distributions 
we find moderate radiative corrections. 
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1 Introduction 



The hard-exclusive photo electroproduction off a nucleoli, among a class of hard exclusive but 
inelastic processes, is generally considered as the theoretic cleanest process to gain access to the 
so-called generalized parton distributions (GPDs) PIE]- These non-perturbative distributions 
contain manifold information that cannot be gained from measurements of deep (semi-)inelastic 
or elastic and other exclusive processes [3]. Based on a partonic interpretation, GPDs allow 
us to reveal the internal structure of the probed hadron, especially of the nucleon, from a new 
perspective, for comprehensive reviews see Ref. j3]. 

We recall that the process in question has two interfering subprocesses, namely, the Bethe- 
Heitler bremsstrahlung and deeply virtual Compton scattering (DVCS). The former is known in 
terms of the electromagnetic form factors and so the latter can be measured in several asymme- 
tries, see e.g., Ref. [2], that appear in leading order of the expansion with respect to the inverse 
photon virtuality, i.e., single and double spin and charge asymmetries jHJ El E]. Relying on 
the validity of the operator product expansion (OPE) of two electromagnetic currents, the DVCS 
amplitude factorizes to leading power, i.e., twist-two, in short- and long-distance physics, where 
the former and latter are incorporated in (resummed) Wilson-coefficients and GPDs, respectively 
PP. Equivalently, in the partonic framework, the DVCS amplitude to leading power accuracy is 
represented as convolution of hard-scattering amplitudes, systematically calculable in powers of 
the strong coupling, and GPDs. The factorization of collinear singularities was shown to all orders 
in perturbation theory |lUj . 

Let us remind that for experimental accessible scales, the perturbative description of elastic 
exclusive processes is controversially debated jllj . With respect to this fact one might wonder 
whether the perturbative GPD framework is justified. For the GPD phenomenology it is essential 
to clarify this issue; unfortunately, this is an interlaced task. The relative size of perturbative and 
non-perturbative, i.e., power suppressed contributions, might serve as criteria for the reliability 
of the perturbative framework. The DVCS hard-scattering amplitudes (Wilson-coefficients) were 
perturbatively evaluated up to next-to-leading order (NLO) in the strong coupling [T2"l I13j. It 
turned out that these radiative corrections can be of the order of 30% to 50% for fixed target 
kinematics [T3]. Even much larger corrections have been reported for the kinematics of collider 
experiments [T^]. Mainly, but not only, this is due the appearances of gluonic GPDs at NLO, see 
also discussion in Sect. 6.2.3 of Ref. The problem arises whether these sizeable corrections 
are related to the appearance of gluons, are induced by an awkward choice of GPD ansatze, or 
indicate that at experimental accessible scales the perturbative regime has not be fully set in. 

To get a deeper insight into this issue, we study in this letter perturbative corrections to the 
DVCS process beyond NLO accuracy within a framework that avoids a cumbersome diagrammatic 
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calculation. Based on the conformally covariant OPE (COPE), the perturbative corrections up 
to next-to-next-to-leading order (NNLO) can be economically evaluated fOJE], see also review 
|17j . Thereby we restrict ourselves to the flavor non-singlet part of the helicity conserved twist- 
two amplitude in the parity even sector, which is the dominant contribution in several DVCS 
observables. In Sect. |21we outline the evaluation of this amplitude by means of the COPE [To! 

This yields a divergent series of conformal GPD moments, which is resummed within a 
Mellin-Barnes integral. In Sect. El we numerically study the radiative corrections, mainly due to 
the Wilson-coefficients, at NLO and NNLO. Especially, we explore the numerical differences of 
the NLO corrections in the modified minimal subtraction (MS) and conformal subtraction (CS) 
schemes. We then provide estimates to NNLO accuracy. Finally, in Sect. |U we conclude. 

2 General formalism 

Let us first recall the standard perturbative QCD framework. Usually, one employs the MS 
scheme, which is based on dimensional regularization and the removal of the poles with respect to 
the dimensional parameter (4 — n)/2. This scheme is used twofold: (i) to define the renormalized 
strong coupling and (ii) for the factorization of collinear singularities or if one wishes for the 
renormalization of (leading twist) composite operators 1 , which are labelled by the good quantum 
number spin. In exclusive reactions such operators with given spin are also be build within total 
derivatives. In general operators with the same spin will mix under renormalization. To get rid of 
this mixing phenomenon, one changes at some stage of the full calculation in an explicit or implicit 
manner to a basis of multiplicatively renormalizable operators, i.e., one diagonalizes the evolution 
equation ^H]- For vanishing /3(a s ) function this transformation leads to conformal operators 
|2Uj . Such operators transform covariantly under the so-called collinear conformal transformation 
SL(2, M) and are members of infinite dimensional conformal multiplets (towers) that are labelled 
by the conformal spin. With other words one takes explicitly or implicitly advantages from the 
underlying conformal symmetry of the classical QCD Lagrangian. If the trace anomaly of the 
energy-momentum tensor, proportional to /3(a s ), is absent, conformal symmetry is present in 
perturbative QCD and its predictive power can be employed at any order. 

In the following we evaluate the scattering amplitude for the DVCS process, which is expressed 
in terms of the time ordered product of two electromagnetic currents 

T tMV (q,P 1 ,P 2 ) = ^J dxe i ^(P 2 ,S 2 \Tj l ,(x/2)j u (-x/2)\P l ,S 1 ), (1) 

1 To leading twist approximation there is an one-to-one cross talk between the ultra-violet and collinear singu- 
larities in the partonic matrix elements of these operators. 



2 



where q = (qi + <?2)/2 (/i and q2 refers to the outgoing real photon). The incoming photon has a 
large virtuality q\ = — Q 2 and requiring that in the limit — q 2 = Q 2 — > oo the scaling variables 

t = -^-, V = -^ L , P = P 1 + P 2 , A = P 2 -P 1 . (2) 
P ■ q P ■ q 

and the momentum transfer squared A 2 are fixed, the dominant contributions arise from the 

light-cone singularities of the time ordered product. In this generalized Bjorken limit one can 

now employ the OPE to evaluate the hadronic tensor in terms of the leading twist-two operators, 

where for DVCS kinematics 77 ~ £ and Q 2 ~ Q 2 /2. 

Before we outline this step let us introduce a parameterization of the hadronic tensor 

q a V a q p A a 
T^viq, P, A) = -g^y— - ii^—— + . . . . (3) 

To ensure current conservation, the metric and Levi-Civita tensors are contracted here with 
projection operators, for explicit definitions of g^ v and e^ vpa see, e.g., Ref. [7]. The ellipsis indicates 
terms that are finally power suppressed in the DVCS amplitude or a determined by the gluon 
transversity GPD, which is not considered here. We note that in the forward limit A — > the 
first and second term on the r.h.s. are expressed by the deep inelastic structure functions Fi and 
gi. In the parity even sector the vector 

v° = u(p 2 , s 2 ) [nr + £^r) u(Pi, s 1 ) + ---, (4) 

is decomposed in Compton form factors (CFFs) 7i and S, similar for the axial-vector A a in terms 
of 7i and S, where again higher twist contributions are neglected. 

Now we are in the position to employ the OPE. Let us first suppose that the trace anomaly is 
absent and so conformal symmetry can be employed. Formally, this can be achieved by assuming 
that a hypothetical fixed point exist, i.e., /3(a*) = 0. We can then use the COPE, which tells us 
how the total derivatives are arranged. Moreover, it can be considered as a partial wave expansion 
with respect to the conformal spin j + 2. For the time-ordered product of two electromagnetic 
currents it reads in the flavor non-singlet sector as |2~Ti ITo] 
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7T 2 (-x 2 + zO) 2 1 K 7J 2^/ 2 r(j + l + 7i /2) 

-2tx.) j+1 [ du [u(l - u)] j+1 ^ 2 Of s (uxJ) + • • • , (5) 
Jo 

where x_ = n ■ x with n 2 = is the projection of the vector x on the light cone. The anomalous 
dimensions of the multiplicatively renormalizable operators 0^ s are denoted as 7^(0;*). The 
advantage of the COPE is that the Wilson-coefficients (electrical charge factors will be omitted) 

c J { x/i)-{ f^x) r(3/2)r(3 + J + 7 J /2) CjK) ' (6) 
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are up to the normalization Cj(a*) known. The conformal operators 
,ns,_ ,_r(3/2)r(l + ;),.. ^j, , _3 /2 /S 




= 2ir(3/2 + j) (*<W(^-h+ ^ J > 

are defined with a non-standard normalization. Here D+=D+ — D+ and d + =D+ + £>+ are the 
covariant and total derivatives, contracted with the light-like vector n, i.e., n 2 = and n ■ n = 1, 
C^ 2 are the Gegenbauer polynomials of order j with index 3/2. The normalization is chosen so 
that in the forward limit A — > the reduced matrix elements of the conformal operators 

-^ TI {P 2l S 2 \O j (0)\Pi^ 1 ) = ^(P,,^)^-^^ (8) 

coincide with the Mellin moments of parton densities. Especially, for odd values of j the forward 
limit of Hj is the sum of the quark and anti-quark density moments, i.e., lim^ Hj = qj + q~j. 
Moreover, in this limit the ^-integration in Eq. (jSJ) can be performed and leads to the well- 
known OPE that is used for the deep inelastic scattering structure function Fx. Hence, Cj(a*) are 
identified as the Wilson-coefficients of Fx, known to NNLO order. 

Employing the COPE JHJ), we can now achieve the factorization of the CFFs (J3J) in a straight- 
forward manner. Plugging Eq. ([""J into the definition of the hadronic tensor (Q), performing Fourier 
transform, and form factor decomposition ("""J, we arrive for the DVCS kinematics rj = £ at [II"] 



(9) 



Here Hj and Ef , cf. Eq. (JHJ), are given by the conformal moments of the corresponding GPDs 




Note that these expectation values are measurable on the lattice ["""""] . 

Unfortunately, the series is divergent for |£| < 1 and must be resummed, e.g., by means of 
the Sommerfeld- Watson transformation Alternatively, using the trick that 77 is not equated to 
£ allows to deal with Eq. J""} in the analogous manner as it is known from deep inelastic scattering. 
Namely, a dispersion relation allows to express the Mellin moments of the imaginary part by the 
partial waves, appearing in the COPE Q, and thus the inverse Mellin transform provides the 
imaginary part of the CFFs. The real part can be restored from the dispersion relation, too, and 
finally equating rj = £ leads to the same representation as presented in Ref. [2*3*] : 

|^ S S } (£, A 2 ,/, 2 ) = j l [ + 2 d ^~ 3 ' 1 [* + tan (y)] Cf s (Q/n) I I (£, A 2 , fi 2 ) , (11) 



where all singularities of the conformal GPD moments and conformal moments lie on the l.h.s. 
of the integration path (c < 1). Note, however, that the analytic continuation of the Wilson- 
coefficients (JBJ) with respect to the conformal spin, analogously done as in deep inelastic scattering, 
leads essentially to an exponential 2 J growing with increasing j. For £ > 1, this must be weighed 
down by a suppression factor that comes from the conformal moments (fTI]|) . This is a rather 
nontrivial requirement on their analytic continuation. It should be done in such a way that the 
integral (jll)) remains unchanged if the integration contour is closed by an infinite arc, surrounding 
the first and forth quadrant. The residue theorem states then that this Mellin-Barnes integral is 
for £ > 1 equivalent to the series @. We note that the form of the integrand does in fact 
not rely on conformal symmetry and so it can be used in any scheme within the corresponding 
Wilson-coefficients. 

In the MS scheme the conformal symmetry is not manifestly implemented in the OPE, however, 
this failure can be cured by a finite renormalization H2] As explained above, in such a 
CS scheme and for j3 = we can simply borrow the Wilson-coefficients Cj(a s ) and anomalous 
dimensions 7j(« s ) from deep inelastic scattering. The inclusion of (3 proportional terms in Eq. 
(jllj) is, as factorization by itself, conventional. Two possibilities have been discussed in Ref. |18j . 
Namely, in the CS scheme, one might add the j3 proportional term from the MS scheme, evaluated 
to NNLO in Ref. [21], while in the CS scheme the running of the coupling is implemented in the 
form of the COPE (j5J) in such a way that the Wilson-coefficients autonomously evolve in the 
considered order [TH|. In this letter we prefer the latter convention and, moreover, will equate the 
factorization and renormalization scales, i.e, [i = fif = /i r . After expansion with respect to a s , we 
write, discarding in the following the superscript NS, the Wilson-coefficients (jUJ) as 
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»f) 3 -W + 3/2) + S 2 (j + 2). (17) 



The analytic continuation of the harmonic sums are defined by Si(z) = d\nT(z + l)/dz + 7^ and 
5*2 ( z ) = ~d 2 m T(z + 1)/ dz 2 + C(2), where 7^ is the Euler constant. The first expansion coefficient 
of (3(g) /g = (a s /47r)/5o + 0(a 2 s ) is /3q = (2/3) rif — 11, where n/ is the number of active quarks. 
The two-loop quantities c? and jj are lengthy and can be obtained from Ref . [23 • The 
evolution of the flavor non-singlet (integer) conformal moments in this CS scheme is governed by 



dp \ Ef 



2tt ~ 7j ' 



(2tt; 



2 



2 (2tt) 3 ^ ' 



fc=0 



(2tt) 3 ' 3 



Ef 



(r/,A 2 ,/i 2 



(^A 2 ,/i 2 



where the mixing matrix A^r is not completely known. 

We remark that the NLO corrections in the MS scheme can be easily evaluated from the con- 
formal moments of the hard-scattering amplitude, e.g., given in Ref. |14j . All integrals, which are 
needed, are given in Appendix C of Ref. ^Hj for integer conformal spin. The analytic continuation 
is straightforward and so in the MS scheme Eq. (J13J) is to replace by 
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5-4^ + 1) 
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(o) 2 



2(j + 1)(j + 2) ( J + l)2( J + 2 )2 
In this scheme also the complete anomalous dimension matrix is known to two-loop accuracy 



(19) 



3 Numerical results 

In this Section we numerically study the radiative corrections to the CFF 

at NLO (P = 1) and NNLO (P = 2) accuracy. Certainly, the size of the radiative corrections 
depend on the GPD distribution itself. As a GPD model assumption, we suppose that the expan- 
sion in powers of £ 2 induces a systematic expansion of the CFF (J2Uj) . A closer look to this 
issue has been given in Sect. 4 of Ref. j2EJ- The GPD moments are generically parameterized as 

rrNS/ A 2 t 2\ _ pNS/ A 2n B (l ~ tt(A 2 ) + j, (3 + 1) 2 

H 3 (A,^ )-F (A B{1 _ a{A2) ^ + 1) +O(0- (21) 

In the forward limit, i.e., A — ► 0, these moments reduce to the Mellin moments of the unpolarized 
parton density. Hence, the parameters a and (3 characterize the small and large x behavior, 
respectively, i.e., q(x,fio) oc x~ a (l — Xf. For non-singlet parton densities the generic values are 
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Figure 1: The relative radiative corrections, defined in Eq. (|22j) . are plotted versus £ for the 
modulus [(a) and (b)] and phase [(c) and (d)] of the CFF (|2U|) with a = 0.5 [(a) and (c)] and 
a = -0.1 [(b) and (d)]: NNLO in full (solid) and p = (dash-dotted) as well as in NLO for the 
CS (dashed) and MS (dotted) scheme. We always set fi = Q and a s (Q 2 = 2.5GeV 2 )/7r = 0.1. 

a(A 2 = 0) = 1/2 and /3 — 3. In the off-forward kinematics we consider a(A 2 ) = a(0) + a'(0)A 2 
as a linear meson Regge trajectory with a(Q) = 1/2 and a'(0) = 0.9GeV~ 2 . Note also that the 
conformal GPD moment with j — reduces to the elastic form factor F NS (A 2 ). Certainly, a more 
realistic ansatz of GPD moments is formed by a linear combination of the building blocks (|21j). 

Let us first compare the radiative corrections in the MS and CS scheme to NLO accuracy. 
Strictly spoken, at a given input scale Qo there is no difference between both predictions, if the 
non-perturbative quantities are transformed, too, and a consequent expansion in a s is performed. 
Usually, the GPD moments are taken from some non-perturbative (model) calculation or ansatz, 
where the matching with the perturbative prescriptions has its own uncertainties. So let us take 
the same input in both schemes and study the relative changes 



K. 



n 



N p LO 



abs 



K 



arg [H 



/N p LO 



(22) 



* arg arg(ft NP - 1LO ) 

to the modulus and phase of the CFF f)20j) . These quantities do not suffer from large radiative 
corrections as it is artificially the case for the real part of the amplitude 2 . One should bear in 
2 The real part possesses a zero in the valence quark region, which position depends on the approximation. 
Hence, in the vicinity of this zero the relative radiative corrections blow up. 
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mind that these factors are a measure for the necessary reparameterization of the GPD when one 
includes the next order in a given scheme. In Fig. we depict for the typical kinematics in fixed 
target experiments, i.e., 0.05 < £ < 0.3, the K factors to NLO as dashed and dotted lines for the 
CS and MS schemes, respectively. We set \i — Q and independent of the considered approximation 
we choose a s ([j% = 2.5 GeV 2 ) = 0.17T. From the panels (a) and (b) it can been realized that in the 
MS scheme the radiative corrections to the modulus are up to 20% and 30% for a = 0.5 (A 2 = 0) 
and a = —0.1 (A 2 = —0.6 GeV 2 ), respectively. In the CS scheme these radiative corrections are 
reduced by 30%. Note that such a reduction has been observed in a quite different processes, 
namely, for the photon-to-pion transition form factor 3 . The relative radiative corrections to the 
phase is in both cases of about 15% at £ = 0.6 and diminishes with decreasing £. These findings 
qualitatively agree with previous ones in which the Radyushkin ansatz for GPDs was used [T%] . 

We study now the radiative corrections to NNLO accuracy. To simplify their evaluation, we 
take for a fit, given in Ref. [27], rather than the exact expression. For the (3q proportional 
term we have checked that within |c| = \^tej\ < 1/2, see Mellin-Barnes integral (|20j) . the accuracy 
is on the level of one per mill or better. Outside this region, the deviation can be larger and one 
might use Fortran routines j2Hj- For three quark flavors n/ = 3, the same scale settings, and the 
initial value of a s as specified above for NLO, we plot in Fig. ^the ratios (J22J) in the CS scheme for 
P = 2 as solid line. The ^-dependence of the modulus i^-factor, see panels (a) and (b), is rather 
flat and the modulus decreases on the 5% level. The radiative correction to the phase is again 
negligible for smaller values of £ and increases now only to 5% at £ = 0.6. We note that the (3$ 
induced corrections are about twice times larger than the remaining ones and are opposite in sign, 
see dash-dotted line. Remarkably, for j3 Q = we find then an opposite behavior of the K factors 
as in NLO. This arise from a sign alternating series and as for the photon-to-pion transition form 
factor, this might be considered as a reminiscence on Sudakov double logs [29J. 

Let us finally address the modification of the scale dependence due to the higher order correc- 
tions. Note that the analysis about scheme dependence in Ref. [3U] suggests that the discrepancy 
between the CS and MS schemes at NLO accuracy are mainly induced by the Wilson-coefficients 
while the evolution yields minor differences. So we only consider here the CS scheme and analo- 
gously as in Eq. (|22[). we quantify the relative changes of dTi/dln Q 2 by the ratios 

Both the numerator and denominator in the latter i^-factor are now defined in the interval [0, 2ir] 

3 This process is also evaluated within the OPE and the reduction of NLO corrections has a common origin. 
Namely, in the CS, compared to the MS, scheme the first few Wilson-coefficients are smaller, while the stronger 
logarithmic growing with increasing j is anyway suppressed by the non-perturbative input, see Eqs. (|13|) and l|f 9|) . 
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Figure 2: The relative change of scale dependence, cf. Eq. ([23]). in the CS scheme at NLO (dashed, 
dotted) and NNLO (solid, dash-dotted) versus £ is depicted for the modulus (a) and phase (b) of 
the CFF (1201) with a = 0.5 (dashed, solid) and a = -0.1 (dotted, dash-dotted) and Q 2 = 4GeV 2 . 
We set n = Q, a s (fi 2 = 2.5 GeV 2 )/vr = 0.1 and took the input (jH} at the scale Q 2 = 0.5 GeV 2 . 

and so the appearance of a zero in the denominator is avoided, cf. footnote El We take the same 
scale setting and initial condition for the (exact) evolution of a s (Q) as above. The conformal 
moments ([21)1 are evolved in the CS scheme. The input scale is now /z 2 , = 0.5 GeV 2 , which is typical 
for the matching of perturbative and non-perturbative QCD. The non-leading logs in the solution 
of the evolution equation (JT5J) are expanded with respect to a s and are consistently combined with 
the Wilson-coefficients (J 12)) in the considered order. Here the forward anomalous dimensions are 
taken from Ref. [31] and the unknown NNLO mixing term A^ fc s in Eq. ()18|) is neglected. This 
mixing can be suppressed at the input scale within an appropriate initial condition and so we 
expect only a minor numerical effect; see also Ref. Ref. [201 • The dashed and dotted lines in Fig- EI 
show that in NLO the scale dependence changes of about 30% to 50% for the modulus and up to 
20% for the phase. The latter iif-factor becomes close to one for smaller values of £. The NNLO 
corrections, compared to the full NLO result, are milder. We observe (solid and dash-dotted line) 
a 20% or less and a maximal 10% correction for the modulus and phase, respectively. 



4 Conclusions 

In this letter we have taken the first steps in the full investigation of radiative NNLO corrections 
to the DVCS scattering amplitude. Thereby, we employed the COPE, which allows an economical 
treatment of perturbative corrections, and restricted ourselves to the CFF 7i in the flavor non- 
singlet sector. This study can be immediately adopted for £. The extension to the axial vector 
case, i.e., to TC and S, is straightforward, however, here the anomalous dimensions to three loop 
order cannot be borrowed from the polarized DIS results. The conformal approach can be also 
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straightforwardly extended to the flavor singlet sector, which will be presented somewhere else. 

We relied on GPDs for which the expansion of the conformal moments in powers of £ 2 is 
meaningful and used in our analysis only the leading term. Whether this assumption is too 
restrictive is an open problem, which might be resolved by lattice calculations [22] • At least there 
is a warning from perturbative QCD. Namely, in the MS scheme ^-suppressed contributions in 
the Wilson-coefficients are resummed and lead to larger radiative corrections than in the CS one. 
Let us summarize our numerical findings for fixed target kinematics, i.e., 0.05 < £ < 0.3: 

• Compared to the MS scheme, in the CS one the radiative NLO corrections to the modulus 
of the CFF are reduced by 30% or so, while for the phase no significant differences appear. 

• The relative NNLO corrections compared to the NLO ones are of the 5% level or below, 
where the /3o proportional ones dominate and determine the sign. 

• The change of the scale dependence due to the NLO corrections is rather large and is for 
the modulus (phase) of the CFF of about 30% to 40% (15% or smaller). Comparing NNLO 
to NLO, the two-loop corrections are reduced to less than 20% (10%) . 

Let us add that the Wilson-coefficients in the vector and axial-vector channel possess a similar 
j dependence. So we might expect that for H and £ (axial-vector case) comparable radiative 
corrections appear. We conclude that our findings support the perturbative description of the 
DVCS process at scales which are accessible in the kinematics of present fixed target experiments. 
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